DUDLEY K170X LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
m .q TTV - ( . vv ^/n.iW'JW 93943-B002 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

HEAT TRANSFER FROM AN ARRAY OF 
HEATED CYLINDRICAL ELEMENTS ON AN 
ADIABATIC CHANNEL WALL 

by 

James D. Piatt, Jr. 

September 1986 

Thesis Adivsor: Allan D. Kraus 

Approved for public release; distribution is unlimited. 



T232219 




UNCLASSIFIED 

SE CuPtir v Classification op This pagT 



REPORT DOCUMENTATION PAGE 



la REPORT SECURITY CLASSIFICATION 



UNCLASSIFIED 



lb RESTRICTIVE MARKINGS 



2a SECURITY CLASSIFICATION AUTHORITY 



2b DECLASSIFICATION /DOWNGRADING SCHEDULE 



3 DISTRIBUTION / AVAILABILITY OF REPORT 

Approved for public release; 
distribution is unlimited 



4 PERFORMING ORGANIZATION REPORT NUM8ER(S) 



5 MONITORING ORGANIZATION REPORT NUMBER(S) 



6a NAME OF PERFORMING ORGANIZATION 

Naval Postgraduate School 



6b OFFICE SYMBOL 
(If applicable) 

Code 69 



7a NAME OF MONITORING ORGANIZATION 

Naval Postgraduate School 



6c ADDRESS (C/fy, State, and ZIP Code) 



Monterey, California 93943-5000 



7b. ADDRESS (City, State, and ZIP Code) 



Monterey, California 93943-5000 



8a NAME OF FUNDING / SPONSORING 
ORGANIZATION 



8b OFFICE SYMBOL 
(If applicable) 



9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 



8c ADDRESS (City, State, and ZIP Code) 



10 SOURCE OF FUNDING NUMBERS 



PROGRAM 


PROJECT 


TASK 


WORK UNIT 


ELEMENT NO 


NO 


NO 


ACCESSION NO 



11 TITLE (Include Security Classification) 

HEAT TRANSFER FROM AN ARRAY OF HEATED CYLINDRICAL ELEMENTS ON AN ADIABATIC CHANNEL WALL 



I 12 PERSONAL AUTHOR(S) 



Ha TYPE OP REPORT 


13b TIME COVERED 


14 DATE Of REPORT (Year, Month, Day) 


15 PAGE COUNT 


Master's Thesis 


FROM TO 


1986 September 


58 


16 SUPPLEMENTARY NOTATION 









17 



COSATI CODES 



FIELD 



GROUP 



SUB-GROUP 



18 SUBJECT TERMS ( Continue on reverse if necessary and identify by block number) 

Heated cylindrical array, convective heat-transfer from 
cylinders, circuit board cooling, electronic cooling. 



9 ABSTRACT { Continue on reverse it necessary and identify by block number) 

This thesis describes an experimental study of the fluid dynamics and heat transfer 
characteristics of Laminar air flow across cylindrical elements mounted on one wall of 
a vertical adiabatic channel. Various combinations of approach velocity and channel 
widths were employed and the variations of the row by row heat transfer coefficients 
between elements of the array and the air stream were determined. 



20 O S TRi3UTiON / AVAILABILITY Of ABSTRACT 

fit] UNCLASSlflEOAJNUMITEO □ SAME AS RPT □ OTIC USERS 


21 ABSTRACT SECURITY CLASSIFICATION 

UNCLASSIFIED 


22a NAME OF RESPONSIBLE INDIVIDUAL 

Prof. Allan D. Kraus 




i i; tom3i<s BOL 



DD FORM 1473, 84 MAR 



83 APR edition may be used until exhausted 
All other editions are obsolete 

i 



SECURITY CLASSIFICATION OF T HI$ PAGE 

UNCLASSIFIED 



Approved for public release; distribution is unlimited 



Heat Transfer From An Array of Heated 
Cylindrical Elements On An Adiabatic Channel Wall 



by 



James D . Piatt, Jr. 

Lieutenant, United States Navy 
B.A., State University of New York at Geneseo, 1980 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN MECHANICAL ENGINEERING 



from the 



NAVAL POSTGRADUATE SCHOOL 
September 1986 



ABSTRACT 



This thesis describes an experimental study of the fluid 
dynamics and heat transfer characteristics of Laminar air flow 
across cylindrical elements mounted on one wall of a vertical' 
adiabatic channel. Various combinations of approach velocity and 
channel width were employed and the variations of the row by row 
heat transfer coefficients between elements of the array and the 
air stream were determined. 



3 



TABLE OF CONTENTS 



I. INTRODUCTION 11 

A. THEORY OF VERTICAL CHANNELS FORMED BY CIRCUIT 

BOARDS 11 

B. THEORY OF SURFACE ROUGHNESS GEOMETRIES 15 

C. THEORY OF HEAT TRANSFER FROM ROUGH SURFACES 15 

II. EQUIPMENT DESIGN AND INSTRUMENTATION 17 

A. WIND TUNNEL 17 

B. TEST SURFACE 19 

C. CYLINDRICAL ELEMENTS 21 

D. INSTRUMENTATION 25 

1. Temperature Measurements 25 

2. Velocity Measurements 25 

3. Power Measurements and Control 26 

4. Data Acquisition 2 6 

E. EXPERIMENTAL PROCEDURE 26 

III. DATA REDUCTION 29 

A. DATA TAKEN 29 

B. SOUGHT AFTER PARAMETERS 29 

1. Dimensionless Temperature Ratio 29 

2. Heat- transfer Coefficient 30 

3. Reynolds Number 31 

4 . Nussel t Number 31 

C. TYPICAL SET OF DATA 32 

D. CALCULATIONS 32 



4 



1. Dimensionless Temperature Ratio 32 

2. Heat-Transfer Coefficient 33 

3. Reynolds Number 34 

4. Nusselt Number 34 

IV. RESULTS 36 

A. THERMAL WAKE 36 

B. HEAT-TRANSFER COEFFICIENT 37 

V. CONCLUSIONS AND RECOMMENDATIONS 47 

APPENDIX A: THERMOCOUPLE CALIBRATION 49 

APPENDIX B: TABULATED DATA 53 

LIST OF REFERENCES 55 

INITIAL DISTRIBUTION LIST 57 



5 



LIST OF TABLES 



I. DATA SET 32 

II. PARAMETERS TESTED 36 



6 



LIST OF FIGURES 



1.1 Modern Circuit Board 12 

1.2 Modular Arrangement of Circuit Board 13 

2.1 Photograph of Wind Tunnel 18 

2.2 Approach Velocity vs. Boundary Layer Thickness for 

Various Channel Lengths 20 

2.3 Sketch of Test Surface Configuration 22 

2.4 Photograph of Cylinders 23 

2.5 Photograph of Cylinders Arranged in an Array on the Test 

Surface 24 

2.6 Photograph of Autodata Nine 27 

4.1 Normalized Dimensionless Temperature vs. Row Position 

for a Channel Heiqht/Cyl i nder Height Ratio of 1.0 With 
Heated Element in the First Row 38 

4.2 Normalized Dimensionless Temperature vs. Row Position 

for a Channel Height/Cylinder Height Ratio of 4.6 With 
Heated Element in the First Row 39 

4.3 Normalized Dimensionless Temperature vs. Row Position 

for a Channel Height/Cylinder Height Ratio of 2.3 and 
Heated Element in the Second Row 40 

4.4 Heat- transfer Coefficient vs. Channel Height/Cylinder 

Height Ratio With No Rows Upstream 41 

4.5 Heat- transfer Coefficient vs. Channel Height/Cylinder 

Height Ratio With One Row Upstream 43 

4.6 Heat- transfer Coefficient vs. Number of Upstream Rows 

for a Channel Height/Cylinder Height Ratio of 1.0 44 

4.7 Heat- transfer Coefficient vs. Number of Upstream Rows 

for a Channel Height/Cylinder Height Ratio of 4.6 45 

4.8 Nusselt Number vs. Reynold Number for all Three Channel 

Height to Cylinder Height Ratios 46 



7 



NOMENCLATURE 



A 


Surface area of cylinder exposed to free stream air, m 2 . 


b 


y intercept. 


C 


Multiplication constant, dimensionless. 


d 


Diameter of cylindrical element, m. 


e + 


Roughness Reynolds number based on e, a ch aracter i s ti c 
height of the roughness, dimensionless. 


g + 


Heat transfer similarity function, dimensionless. 


h 


Local heat- transfer coefficient, W/m 2 -°C. 


H 


Distance between channel splitter plates m. 


I 


Cur rent d rawn , A . 


k 


Thermal conductivity of free stream air, W/m-°C. 


L 


Height of cylindrical element, m. 


m 


SI ope of a line. 


N 


Row position number in array. 


Nu h 


Nusselt number, hH/k. 


Pr 


Prandtl number, C p p/k, dimensionless. 


Q l np ut 


Power input from electrical power supply, W. 


Ql os ses 


Power lost from cylindrical elements, W. 


Re 

H 


Reynolds number, HU<»/v > dimensionless. 


Ts »o 


Surface temperature of the cylinder heated by the 
resistor, °C. 


T S .1 


Surface temperature of the ith cylinder heated by the 
thermal wake, °C. 
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T Freestream air temperature, °C. 

oo 

U Freestream air velocity, m/s. 

00 

V Power supply voltage, V. 

Greek Symbols 

6 Dimensionless temperature of a cylinder in the same . 

column as the cylinder heated by a resistor. 

v Kinematic viscosity, m^/s. 
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I. INTRODUCTION 



The removal of heat from electronic equipment to assure 
reliable operation has been a problem plaguing the electronic 
industry for many years. As a result, many researchers have - 
attempted to develop specific correlations to enable designers to 
attempt to predict possible solutions to this problem. 

This research concerns air flowing across cylindrical 
elements mounted in an array on one wall of a vertical adiabatic 
channel. This configuration is much the same as an array of 
transistors in the T05 configuration arranged on stacks of 
parallel circuit boards in various types of electronic equipment. 

A. THEORY OF VERTICAL CHANNELS FORMED BY CIRCUIT BOARDS 

Many electronic systems consist partially of arrays of 
printed circuit boards supporting a fairly dense array of 
electronic elements on one side as shown in Figure 1.1. These 
circuit boards can be arranged so that the side with the elements 
(front side) faces the relatively smooth backside of another. 

When a series of circuit boards are arranged in such a fashion 
(see Figure 1.2), long vertical channels in which cooling air can 
be supplied by forced convection can result. During operation, 
the elements such as flatpacks, dips resistors and transistors 
generate heat, the majority of which is dissipated to the 
surroundings with little transmitted back through the circuit 
board. As a result, it appears as if all of the heat is 
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Figure 1.1 Modern Circuit Board 
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Figure 1.2 



Modular Arrangement of Circuit Board 
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transmitted by one side of the channel or one side of each 
circuit board to the air stream or air layer between the boards. 

Complexity is added to the system when elements serve 
different purposes and, as a result, possess different sizes, 
shapes and heat dissipating character!' sties . 

The process of heat removal from electronic systems creates a 
very comlex design problem governed by such factors as circuit 
board/element arrangement, spacing, heating characteri sti cs and 
operating environment. Moreover, the ultimate goal of any 
thermal analysis and the prediction of operating temperatures is 
the assessment of reliability. Reliability is the probability 
that an electronic compo nent/ sys tern will operate properly without 
failure for its intended lifetime under its environmental 
operating conditions. [Ref. 1] 

For many years, there have been numerous attempts to derive 
techniques to aid designers in this field. Unfortunately, the 
methods developed are often too complex, inadequate or restricted 
to specific conditions. 

Early work with natural convection between parallel arrays of 
vertical circuit boards was performed by Elenbaas [Ref. 2], 

Bodoia and Osterle [Ref. 3] and Engel and Muller [Ref. 4]. In 
1972, Aung et al . [Ref. 5] studied the effects of circuit board 
supports and elements in the field of flow for wide channel 
spacings. Aung went on to investigate uniform heat flux on two 
channel walls in which he subsequently studied the effects of 
baffles and circuit board staggering. [Ref. 6] There is 
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considerable literature describing the efforts that have been 
made to attempt to model circuit channel cooling by natural 
convection and forced convection. 

For example, Kelleher utilized a thermal network approach to 
derive a more representative model of forced convection for o.ne 
specific geometry with equal heat releases [Ref. 7]. Other 
experimental or semi -expe rime ntal approaches were performed by 
Tamura [Ref. 8], Finch [Ref. 9], Laermer [Ref. 10], Marto [Ref. 
11] and Parkes and Preston [Ref. 12]. Most of this work is 
restricted to specific conditions. 

B. THEORY OF SURFACE ROUGHNESS GEOMETRIES 

Surface roughnesses have been classified by several different 
techniques of which the most commonly used are the k and d type 
roughness classifications consisting of 2 and 3 dimension classes 
which are further divided into artificial and commercial 
roughness subclasses. Although many researchers allude to these 
classifications, only a limited number actually use this 
terminology in agreement. It is further noted that in the area 
of heat- transfer, the k- and d- type roughnesses seem irrelevant 
because the same heat transfer roughness parameter (g + ) 
correlation is proposed for all types of roughnesses. 

C. THEORY OF HEAT TRANSFER FROM ROUGH SURFACES 

Earlier researchers in this field seem to agree that the 
heat- transfer roughness parameter (g + ), commonly referred to as 
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the inverse of the sublayer Stanton number, takes the functional 
f o rm : 

P Q 

9 + = 1 /S t s 1 ay e r = C(e ) (Pn) 

where e + is the roughness Reynolds number based on an arbitrary 
characteri Stic height of the roughness (e). The different va.lues 
of p and q fall into two categories: 

1. The assumption of similarity near the wall which results 
in value of p = .2 and q = .44 or for the eddy-di f f u s i v i ty 
models p = .25 and q = .66 [Refs. 13 and 14]. 

2. The assumption of a breakdown in close to the wall 
similarity results in values of p = .45 and q = .57 or 
p = .5 and q = .66 for the eddy diffusivity models 
[Refs. 15 and 14] . 

With regard to this study, the foregoing is only of academic 
interest because the geometry here is defined as an in-line array 
of cylindrical elements. 
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I I . EQUIPMENT DESIGN AND INSTRUMENTATION 



The equipment utilized in this research was either 
manufactured by the Engineering Machine Shop or supplied by the 
Naval Engineering Department of the Naval Postgraduate School., 
Monterey, California. 

A. WIND TUNNEL 

The wind tunnel (Figure 2.1) is constructed of 1/8" 
plexiglass with the exception of two 1/4" access plates and 
various reinforcement members mounted on the exterior so as to 
ensure smooth air flow on the tunnel interior. 

The test section of the tunnel has a 12" by 4" cross section 
and a height of 18" with access covers on both faces to 
facilitate instrumentation and support for the test surfaces. 

The tunnel intake is belled with an 18" radius to produce an 
orifice that has a 40" by 12" cross section. An adjustable leg 
is attached at each corner to serve as a support for the tunnel, 
to level the tunnel and to provide an 8" clearance off the floor. 

The upper end of the tunnel is reduced to a section of 2" PVC 
ducting which is connected to a 550 cfm turbo fan. Air flow 
through the tunnel is controlled by a series of valves connected 
in 1 i ne with the fan. 

A 1/8" plexiglass splitter plate was mounted 1/2" off of the 
access plate facing the test surface in order to reduce the 
effects of the boundary layer formed in the tunnel. This 
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Figure 2.1 



Photograph of Wind Tunnel. 



18 






distance was pre-de termi ned through the use of computer software 
[Ref. 16] and the results are displayed in Figure 2.2. A wooden 
frame box (4'x4'x 8') was constructed around the wind tunnel to 
elliminate the effects of incident radiation on the cylinders 
from the laboratory environment and the box also served to reduce 
any transient air currents created by movement in the room (i.e., 
doors opening and closing). 

B. TEST SURFACE 

The design of the test surface used in this research was 
derived from the previous work of Arvizu and Moffat [Ref. 17] and 
Ortega [Ref. 18]. It was constructed of a 12" x 12.5" x 1/8" 
sheet of plexiglass (for support) with a 12" x 12.5" x 3/8" sheet 
of balsa wood (thermal boundary) cemented to it. The leading 
edge of this (balsa wood-plexiglass) splitter plate was cut at a 
35° angle to form a sharp leading edge to reduce the effects of 
leading edge separation. The surface (balsa wood) of the 
splitter plate was then sanded and painted three times with white 
enamel paint and sanded a fourth time to ensure a smooth test 
surface. Small (1/16" diameter) holes were drilled through the 
plate to permit the installation of resistor leads and 
thermocouples into the cylindrical elements mounted on the balsa 
wood. Three mounting screws were attached to the plexiglass side 
of the splitter plate and to the access plate to facilitate the 
positioning of the test surface. A tail fin mounted to the tail 
end of the splitter plate was used to regulate the air flow on 



19 



o 

o 




(s/w) iioij aiv 



00 

00 

<D 

c 

u 



S- 

<D 

CO 



>> 

s- 

CO • 
TD to 

c _c 

3 -M 
O CD 
CQ C 
CD 



00 

> i — 

CD 
>) C 
4 -> C 
•r- ra 

o jc 

o o 



CD oo 




U S- 
n3 co 
O > 

s- 

CL S- 
CL O 
< < 4 - 



CNJ 

(XI 



a> 

s- 

3 

CD 



20 



both side of the splitter plate. A fourth adjusting screw was 
attached to this fin and through the access plate 
(Figure 2.3). 

C. CYLINDRICAL ELEMENTS 

The cylindrical elements used in this research were 
manufactured from 3/4" round aluminum stock and cut to a height 
of 1.0 cm. A 7/16" diameter hole was cut to a depth of 1/4" in 
the center of the element and two additional holes (1/8" 
diameter, 1/4 deep) were cut opposite each other to permit the 
installation of the 56 ohm, 1/2 watt resistors. Another hole 
(1/32" diameter, 3/8" deep) was drilled to position the 
thermocouple (Figure 2.4). 

The resistors were epoxied into forty two of the cylinders 
and six more were cut in half so that an array of 7x6 cylinders 
could be mounted evenly on the splitter plate and the 12 half 
cylinders could be arranged along the edges of the plates (Figure 
2.5) . 

The resistor (56 ohm, 1/2 watt), prior to insertion into the 
cylinders, had the leads cut and smaller leads (2841/1, #30 solid 
spc.) attached so as to reduce heat losses by conduction through 
them. 

The cylinders were mounted on the test surface in an in-line 
array with a separation of one cylinder diameter between 
cylinders. This was done both vertically and horizontally. 
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Flow Adjustment 




Figure 2.3 Sketch of Test Surface Configuration. 
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Top Row: Assorted Views of Cylinders 

Bottom Row: T05 Transistor and a nickel 



Figure 2.4 



Photograph of Cylinders. 
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Figure 2.5 Photograph of Cylinders Arranged in an Array 
on the Test Surface. 
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D. INSTRUMENTATION 



1 . Temperature Measurements 

Type T, cop per-constantan thermocouples were utilized for 
the determination of the temperature of the heated elements and 
the inlet and outlet temperatures of the wind tunnel near the 
test surface. The thermocouples were connected to seven scanner 
boards in an autodata nine and calibrated to an accuracy of 
+0.2°F (Appendix A.) The thermocouples were then routed through 
an access plate in the tunnel through the splitter plate and into 
each of the elements. 

2 . Velocity Measurements 

The velocity measurements were obtained by the use of a 
micro manometer and a pitot-static tube. Changes in the dynamic 
and static pressure were measured by the micro manometer and 
using Barnoulli's equation, velocity was determined. The pitot- 
static tube was placed between the splitter plate and the test 
surface splitter plate and used to determine the approach 
velocity between them. 

A TSI (Model 1050/1050AA) hot wire anemometer was used to 
determine the velocity profile prior to entry into the test area. 
The hot wire anemometer was calibrated with the micro manometer 
and used only to ensure that the approach velocity was uniform. 

Accuracy of the mi rco manometer was determined to be 
+_.0005 inches of water and resulted in an error of + 0.2 ft/sec or 
+7cm/s. Because the micro manometer was utilized for the 
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calibration of the hot wire anemometer, more accuracy might 
result if the micro manometer was used to determine the velocity 
between the splitter plates. 

3 . Power Measurements and Control 

A LAMBDA (60 volt, 10 amp) Regulated Power Supply (Model 
LK 345AFM) was used to power to the 56 ohm, 1/2 watt resistors 
mounted in the forty two cylindrical elements. 

A switch board was set up to allow power to be delivered 
to individual elements in either of two columns, the third or the 
fourth columns. It was also set up in such a manner that each or 
all rows could be energized. 

A Hewlett Packard digital multimeter (3466A) was used to 
measure the voltage drop across the resistors and a Fluke 8060A 
multimeter was used to measure the current through the circuit. 
The power delivered is the product of current and voltage. 

4 . Data Acquisition 

The Autodata Nine mentioned above is a self contained 
recording device with up to 100 channels available for both 
temperature and voltage measurements. During this research, 71 
channels were utilized, one for voltage measurement and the rest 
for temperature measurements in degrees Fahrenheit. (Figure 2.6) 

E. EXPERIMENTAL PROCEDURE 

For this research, data was recorded at three different 
channel height to cylinder height ratios (1.0, 2.3 and 4.6) for 
two different approach velocities (1.5 m/s and 3.0 m/s). The 
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Figure 



2.6 Photograph of Autodata 



Nine* 
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data recorded consisted of the voltage supplied to the circuit, 
the current through the circuit and the temperatures of the test 
elements as well as the inlet and exhaust temperatures of the 
airstream. For each channel height and approach velocity, data 
was gathered after a steady state condition was reached, for 
each element of a column energized one at a time and for the 
first four rows and for all six rows energized all at once. 
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III. DATA REDUCTION 



A. DATA TAKEN 

Inlet Free Stream Temperature, T (°F) 

Cylinder Surface temperature, T s , 0 (°F) 
Atmospheric Pressure, P a tm (inches Hg.) 
Pitot-Static Pressure Difference, AP (inches H 2 O) 
Power Supply Voltage, V (volts) 

Current Drawn, I (amps) 



The Dimensionless Tempe ra ture Ra ti 0 , 0, is obtained as a 



function of cylinder location, channel configuration and 
freestream velocity. Here 0 is defined as 



where T s ,i is the surface temperature obtained by the first 
passive cylinder directly after the cylinder heated by a 
resistor. T s ,-j is the surface temperature obtained by individual 
passive cylinders after and in the same column as the cylinder 
heated by the resistor where i is the row after the heated 
cylinder. Note that 




B. SOUGHT AFTER PARAMETERS 



1 . Dimensionless Temperature Ratio 




00 



00 




00 



00 



= 1.0 by definition 
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2. Heat- trans fer Coefficient 



The Heat-transfer Coefficient, h ( W/m 3 °c ) , is found and 
is a function of the cylinder location, channel configuration and 
freestream air velocity. 

The heat transfer coefficient was determined as follows: 

q (W) - Q (W) 

h= input losses 

A ( T s , o - TZ ) 

T s , 0 and T^ were recorded by the autodata nine and 
corrected by the calibration data in Appendix A. 

A, is the surface area of the cylinder is calculated for 
channel /cyl i nder height ratios of 1.0, 2.3 and 4.6. If the 
channel /cyl i nder height ratio is equal to 1.0, A = udL and if the 
c hannel /cy 1 i nde r height ratio is greater than 1.0, A = tt [ d L + 
d 2 /4] . 

The value of Qiosses 1n watts is derived from a parabolic 
temperature dependent approximation utilized by Ortega [Ref. 18] 
in his researach with aluminum cubes. This approximation is 
assumed to be the same for cylinders and is used here. 

Olosses ( W ) = "0,62927 x 10‘ 3 + 0.5059 x 1 0 -2 ( T s , q-T^ ) 

- 0.023695 x lO- 5 (T S , 0 -TJ 2 

Ortega also determined the radiation losses to be an order of 
magnitude smaller than those given above, and radiation losses 
were therefore disregarded. 



30 



Qinput( w ) 1 s the Product of power supply voltage and 
current drawn. 

Qinput = V(volts) x I(amps) 

3 . Reynolds Number 

The Reynolds number is developed as a function of channel 
configuration and free stream air velocity. In this research, 
the Reynolds number was calculated using the freestream air 
velocity, channel height and viscosity as follows: 

Re = channel height x velocity 
H vi sc os i ty 

The value of the viscosity is obtained by a linear 
interpolation of data from tables of viscosity and air 
temperature. 

4 . Mussel t Number 

The Mussel t number is prescribed as a function of channel 
configuration. Here, the Nusselt number was calculated using 
channel height, heat- transfer coefficient and the thermal 
conductivity of the free stream air as follows: 



Mu 



H 



channel height x heat-transfer coefficient 
thermal c o n d u c t i v i ty 



The value of the thermal conductivity was obtained by a 
linear interpolation of data from tables of thermal conductivity 
and air temperature. 



31 



C. TYPICAL SET OF DATA 

Table I presents a typical set of data obtained. These data 
will be used here to illustrate the data reduction procedure. 



p atm = 29.736 inches Hg. 

AP = 0.0051 inches H 2 O 
Channel/ Cylinder Height Ratio = 4.6 
Cylinder Heated - row 1, column 3 
Cylinder Surface Temperature, T s ,q = 206. 1°F 
Voltage Supplied = 14.380 volts. 

Current Drawn = 0.225 amps. 



D. CALCULATIONS 

1 . Dimensionless Temperature Ratio, 

The temperatures for the cylinder surface were recorded 
by the autodata nine and corrected by use of the thermocouple 
calibration data presented in Appendix A. 

T s , (corrected) = m x T s ,-j (recorded) - b 
m Recorded b^ ^s , i (Co rrected ) 

T s>1 = 1.0026 (80.3) - 0.9038 = 79.6 °F 



TABLE I 



DATA SET 




80 . 3°F 
74 . 0°F 
72.2°F 
71 .5°F 
70.8°F 
68 . 6°F 
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T s »2 


= ( 


D. 99994 


(74.0) - 


0.19287 = 73.8 °F 


T S ,3 


= 1.0018 


(72.2) - 


0.05645 = 72.3 °F 


T s .4 


= 1,0018 


( 68 . 6 ) - 


0.5980 = 68.1 °F 


0 


1 

i 


r s >i - T 


00 






i 


r s»l - T 


00 




01 = 


79.6 


- 68.1 


= 11.5 


= 1.0 (by definition) 




T37TT 


- 68.1 


rmr 




6 2 = 


73.8 


- 68.1 


= 5.7 


= 0.4956 




797IT 


~ 68 '. 1 


TT3 




9 3 = 


72.3 


- 68.1 


= 4.2 


= 0.3652 




777T 


- 68.1 


TITS' 




6 4 = 


71.2 


- 68.1 


= 3.1 


= 0.2696 




7TT6 - 


- 68.1 


TITS' 




6 5 = 


70.3 


- 68.1 


= 2.2 


= 0.1913 




79.6 


~627T 


TT7T 




2 . 


Heat-Tansfer 


Coefficient 



h = Q i nput - Ql osses 
A ( T ^ , Q ~ T co ) 

T s , 0 = 1.0025 (206. 1°F) - 0.4745 = 206. 1°F (corrected) 

= 1.0018 ( 68 . 6° F ) - 0.5980 = 68.1°F (corrected) 

T s , 0 - To, = 206.1 - 68.1 = 138. 0°F = 76.7°C 

Channel /cyl i nder height ratio = 4.6 
therefore A = ir[(0.019m) (0.01m) + ( 0 .019m) 2/4] 
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A = 0.00088m 2 

Qlosses = 0-6927 x 10-3 + 0.5059 x 1 0" 2 ( T s , q-T^ ) 

- 0.23695 x 10- 6 (T s , o -TJ 2 

O * 

where T s>0 -Too is in C and Qlosses 1 ' s i n watts. 

T s , 0 -Too — 76.7 C 

Q 1 o s s g s * 0.373 W3 tt s 

• 

Oinput = VxI = l 4 -380 volts x 0.225 amps = 3.27 watts 

h = 3.27 - 0.373 = 42.92 W/m 2 °C 

0700088 (76.7) 

3 . Reynolds Number 
Re H = HU ro / v 

H = 4. 6 ( .01m) = 0.046m 
0 oo — 1 • 5m/s 

= 38.02 x lO- 6 - 109.8 x 10~ 9 (27.0 - T w ) 
where T is in °C 

oo 

T = 20.0°C 

oo 

= 38.02 x 1 0- 6 m 2 / s 

Re^ = .046m x 1.5m/s = 1814.7 

38.02 x 10~ 6 m 2 /s 

4 . Nusselt Number 
Nu h = Hh/k 

H = .046m h = 42.92 W/m 2 °C 

k = 0.0407 - 76.8 x 10“ 6 (27.0-T ) 

OO 
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where T 



oo 



is in 



o 



c 



k = 40.162 x 10- 3 w/m°C 

Nu h = .046m x 42.92 W/m 2 °C = 49.16 
40.162 x IQ" 3 W/m°C 
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IV. RESULTS 



The experimental measurements that were obtained lead to two 
different subject categories: thermal wake effects and 

heat-transfer coefficient. Table II presents the range of the 
parameters covered and a listing of the data is tabulated in 
Appendix C. 



TABLE II 

PARAMETERS TESTED 



Thermal Wake 


Channel Height/Cylinder Height Ratio 
Downstream Row Number 
Vel oci ty (m/s) 


1.0 

0,1, 2, 5 
1.5, 3.0 


4.6 

0,1, 2, 5 
1.5, 3.0 



Heat-Transfer Coefficient 



Channel Height/Cylinder Height 


Rati o 1.0 


2.3 


4.6 


Heated Cylinder Row Number 


1,2, 3, 6 


1,2, 3, 6 


1,2, 3, 6 


Number of Rows Heated 


4,6 


4,6 


4,6 


V e 1 o c i ty ( m / s ) 


1.5, 3.0 


1.5, 3.0 


1.5, 3.0 


A. THERMAL WAKE 

The thermal wake created by 


a single heated cylinder 


in an 



array is influenced by the channel height to cylinder height 
ratio, the mass flow rate and the position of the cylinder in 
the array. 

For a channel height/cylinder height ratio of 1.0, there, is 
no free stream, and the flow through the array undergoes periodic 
accelerations and decelerations. It therefore seems unlikely 
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that a channel height/cylinder height r atio of 4<6 W0U 1d possess 
the same convective means for determination of the downstream 
wake. However, a plot of 0 as a function of downstream row 
number (Figure 4.1 and 4.2) demonstrates a 1 / N dependence for 
both velocities at both channel heights. 

Figure 4.3 shows an even closer resemblance to the 1 / N curve 
for the thermal wake created by a heated cylinder located in a 
row internal to the matrix of cylinders. 

The following is provided for the understanding of the 
graphical curves found in the legends of various graphs in 
Chapter IV : 

Hx . xRxVx . x 

Hx.x Channel Height/Cylinder Height Ratio x.x 

Rx Row that the heated cylinder is in x 

Vx.x Approach Velocity x.x 

For example, H.4.6R1V3.0 denotes data gathered at a channel 
height/cylinder height ratio of 4.6, where the cylinder heated 
is in the first row and the approach velocity is 3.0 m/s. 

B. HEAT-TRANSFER COEFFICIENT 

In Figure 4.4 the heat- transfer coefficient, h, for a single 
heated cylinder in the first row is plotted against channel 
height/cylinder height ratios for approach velocities of 1.5 m/s 
and 3.0 m/s. It can be seen that as the channel height 
increases, the value of h drops off fairly rapidly to an 
asymototic limit at a channel height/cylinder height ratio of 
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V13HI 
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Figure 4.1 Normalized Dimensionless Temperature vs. Row 

Position for a Channel Height/ Cylinder Height Ratio 
of 1.0 With Heated Element in the First Row. 
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Figure 4.3. Normalized Dimensionless Temperature vs. Row 

Position for a Channel Height/Cylinder Height Ratio 
of 2.3 and Heated Element in the Second Row. 
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Figure 4.4 Heat- trans f er Coefficient vs. Channel Height/ 
Cylinder Height Ratio With No Rows Upstream 



4.6. Air flows over the cylinders instead of through them (the 
path of least resistance). 

Figure 4.5 show a plot of h against the channel height to 
cylinder height ratio for an internal cylinder where h drops off 
a little more rapidly. 

Figure 4.6 plots h against array position for a channel 
height/cylinder height ratio of 1.0. Here, h increases with 
position upstream. This correlates to the increased turbulence 
with position downstream. 

Figure 4.7 is a plot of h against array position for a 
channel height to cylinder height ratio of 4.6. For this 
configuration h appears to increase slightly and then drop off 
slightly so that the value of h at the sixth row is less then 
that of the first row. 

Figure 4.8 plots Nusselt number against Reynolds number for 
all three channel height to cylinder height ratios. This 
composition leads to a correlation of the form Nu = C Re0*55 
where C is different for different array densities as noted by 
Arvizu and Moffat [Ref. 17]. 
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Figure 4.5 Heat- transfer Coefficient vs. Channel Height/ 
Cylinder Height Ratio With One Row Upstream. 
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Figure 4.7 Heat- transfer Coefficient vs. Number of Upstream 
Rows for a Channel Height/Cylinder Height Ratio 
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Figure 4.8 Nusselt Number vs. Reynold Number for all Three 
Channel Height to Cylinder Height Ratios. 



V. CONCLUSIONS AND RECOMMENDATIONS 



In this research, for cylindrical elements mounted in a 
vertical channel, the characteristic curves produced closely 
resemble those for a similar configuration of cubical elements 
which were investigated by of Arvizu and Moffat [Ref. 17]. 

Although there is room for more data to assist in making firm 
conclusions about the flow characteristics of air flow around 
cy 1 i nde rs , the 

The following has been observed: 

1. The 1/N characteristic curve produced by the thermal wake 
of a heated cylinder on an adiabatic wall in a channel, 
may allow for fairly accurate prediction of temperature 
increase in downstream elements in an array for various air 
f 1 ow rates . 

2. Channel spacings and airflow rates have very little affect 
on the shape of the characteristic curve produced by the 
thermal wake of a heated cylinder in an array. 

3. The position of a heated cylinder in an array, on an 
adiabatic wall in a channel, has a definite affect on its 
heat transfer characteri sties . Elements requiring more 
cooling should be placed closer to the freestream air inlet. 
This would reduce effects created by the thermal wake of 
other cylinders. 

4. Heat transfer in an array of cylindrical elements is 
mainly influenced by the air flow between the cylinders 
and slightly less by the turbulence in the air flow. 

5. As the channel height to cylinder height increases, the 
amount of heat transfer that occurs decreases to an 
asymototic limit near a channel /cyl i nder height ratio of 
4.6. 

For continued work in this field, it is recommended that, a 
longer test surface with more rows of elements be utilized along 
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All of this 



with higher velocities and various geometric shapes, 
can be accomplished with only a few modifications to the 
equipment at hand. 
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APPENDIX A 



THERMOCOUPLE CALIBRATION 

The thermocouples were all connected to a wooden framed 
switchboard which was connected by 10' wire leads to seven 
scanner boards mounted in the autodata nine data acquisition 
device. This combination was then calibrated as a complete unit 
in the calibration laboratory against five known temperatures: 
66.02°F, 1 06 . 1 4°F , 127.615°F, 151.20°F and 200.42°F 
Thermocouple readings follow: 



1 


66.3 


106.4 


127.9 


151.3 


200.5 


35 

36 


66.3 

66.3 


106.5 

106.5 


127.9 

128.0 


151.2 

151.4 


200.4 

200.5 


2 


66.3 


106.5 


128.0 


151.3 


200.6 


37 


66.3 


106.5 


127.9 


151.3 


200.5 


3 


66.2 


106.5 


128.0 


151.2 


200.4 


38 


66.4 


106.5 


128.0 


151.3 


200.6 


4 


66.3 


106.5 


128.0 


151.3 


200.6 


39 


66.3 


106.4 


127.9 


151.2 


200.4 


5 


66.1 


106.4 


127.9 


151.2 


200.4 


40 


66.6 


106.7 


128.2 


151.6 


200.7 


6 


66.2 


106.4 


127.9 


151.2 


200.4 


41 


66.5 


106.6 


128.1 


151.5 


200.7 


7 


66.1 


106.2 


127.8 


151.1 


200.3 


42 


66.6 


106.8 


128.4 


151.6 


200.9 


8 


66.2 


106.4 


127.8 


151.3 


200.4 


43 


66.6 


106.7 


128.2 


151.5 


200.8 


9 


66.1 


106.2 


127.7 


151.0 


200.2 


44 


66.7 


106.8 


128.4 


151 .6 


200.8 


10 


66.2 


106.4 


127.8 


151.2 


200.4 


45 


66.5 


105.8 


128.2 


151.5 


200.7 


11 


66.8 


106.9 


128.5 


151.8 


201.0 


46 


66.6 


106.8 


128.3 


151.6 


200.8 


12 


66.8 


107.0 


128.5 


151 .8 


200.9 


47 


66.5 


106.7 


128.2 


151.6 


200.7 


13 


66.7 


106.8 


128.3 


151.6 


200.8 


48 


66.7 


106.7 


128.2 


151.6 


200.7 


14 


66.7 


107.0 


128.4 


151.8 


201.0 


49 


66.6 


106.6 


128.0 


151.5 


200.7 


15 


66.6 


106.8 


128.3 


151.6 


200.7 


50 


66.4 


106.5 


128.1 


151 .4 


200.5 


u 


66.7 


106.8 


128.3 


151.7 


200.8 


51 


66.3 


106.4 


127.9 


151.3 


200.6 


17 


66.5 


106.7 


128 .0 


151.4 


200.7 


52 


66.4 


106.6 


128.2 


151.4 


200.7 


18 


66.7 


106.8 


128.2 


151.6 


200.8 


53 


66.3 


106.6 


128.1 


151.4 


200.6 


19 


64.5 


106.6 


128.0 


151.4 


200.7 


54 


66.4 


106.6 


128.2 


151.4 


200.7 


20 


64.2 


106.8 


122.6 


152.4 


199.3 


55 


66.4 


106.6 


128.0 


151.4 


200.6 


21 


66.2 


106.4 


128.0 


151.3 


200.4 


56 


66.4 


106.6 


128.3 


151.5 


200.7 


22 


66.3 


106.5 


128.0 


151.3 


200 .5 


57 


66.4 


106.5 


128.0 


151.4 


200.6 


23 


66.1 


106.3 


127.9 


151.2 


200.4 


58 


66.6 


106.6 


128.2 


151.4 


200.6 


24 


66.2 


106.5 


128.0 


151.3 


200.5 


59 


66.5 


106.6 


128.1 


151.4 


200.6 


25 


66.0 


105.3 


127.8 


151.2 


200.3 


60 


66.2 


106.8 


128.0 


151.3 


200.4 


26 


66.2 


106.5 


127.8 


151.2 


200.4 


61 


66.2 


106.7 


128.0 


151.4 


200.5 


27 


66.0 


106.2 


127.1 


151.0 


200.2 


62 


66.3 


107.1 


128.2 


151.4 


200.6 


28 


66.1 


106.2 


127.7 


151.1 


200.2 


63 


66.3 


106.8 


128.0 


151.3 


200.6 


29 


66.0 


106.2 


127.7 


151.1 


200.3 


64 


66.3 


106.9 


128.2 


151.4 


200.6 


30 


66.3 


106.5 


127.8 


151.2 


200.1 


65 


66.3 


106.8 


128.1 


151.4 


200.6 


31 


66.3 


106.4 


127.8 


151.2 


200.4 


66 


66.3 


106.8 


128.1 


151.4 


200.6 


32 


66.3 


106.5 


128.0 


151.4 


200.6 


67 


66.4 


106.9 


128.3 


151.5 


200.7 


33 


66.2 


106.3 


127.9 


151.4 


200.6 


68 


66.3 


106.8 


128.1 


151.5 


200.7 


34 


66.3 


106.5 


128.0 


151.4 


200.6 
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A least squares fit program was developed to produce 
thermocouple corrections with the following result: 
Tactual = M x Treasured + 8 



THERMOCOUPLE 


i 


M» 


1.0017233 


THERMOCOUPLE 


2 


M» 


1.0012808 


THERMOCOUPLE 


3 


M» 


1 • 0022497 


THERMOCOUPLE 


4 


M* 


1.0012808 


THERMOCOUPLE 


5 


Ma 


1.0013437 


THERMOCOUPLE 


4 


Ha 


1.0019844 


THERMOCOUPLE 


7 


Ma 


1.0017433 


THERMOCOUPLE 


B 


Ma 


1.0017471 


THERMOCOUPLE 


9 


Ma 


1.0024321 


THERMOCOUPLE 


10 


Ma 


1.0019540 


THERMOCOUPLE 


11 


Ma 


1.0017424 


THERMOCOUPLE 


12 


Ma 


1.0024770 


THERMOCOUPLE 


IS 


Ma 


1.0024207 


THERMOCOUPLE 


14 


Ma 


1 . 00131B0 


THERMOCOUPLE 


15 


Ma 


1.0024B44 


THERMOCOUPLE 


14 


Ma 


1.0024144 


THERMOCOUPLE 


17 


Ma 


1.0021400 


THERMOCOUPLE 


IB 


Ma 


1.0025959 


THERMOCOUPLE 


19 


Ma 


1 . 001B959 


THERMOCOUPLE 


20 


Ma 


0.9914704 


THERMOCOUPLE 


21 


Ma 


1.0018024 


THERMOCOUPLE 


22 


Ma 


1.0019808 


THERMOCOUPLE 


23 


Ma 


1.0011015 


THERMOCOUPLE 


24 


Ma 


1 . 00133B0 


THERMOCOUPLE 


25 


Ma 


0.9984798 


THERMOCOUPLE 


24 


Ma 


1.0022001 


THERMOCOUPLE 


27 


Ma 


1 . 001B101 


THERMOCOUPLE 


2B 


Ma 


1.0024242 


THERMOCOUPLE 


29 


Ma 


1.0010777 


THERMOCOUPLE 


30 


Ma 


1.0049429 


THERMOCOUPLE 


31 


Ma 


1.0024007 


THERMOCOUPLE 


32 


M» 


1.0010719 


THERMOCOUPLE 


33 


Ma 


0.9999JB0 


THERMOCOUPLE 


34 


Ma 


1.0010719 



Ba 


•0.4258300 


THERMOCOUPLE 35 


Ba 


•0.4281250 


THERMOCOUPLE 34 


B« 


•0.4744428 


•THERMOCOUPLE 37 


Ba 


-0.4281250 


' THERMOCOUPLE 38 


Ba 


-0.2941913 


THERMOCOUPLE 39 


Ba 


•0 • 399755B 


THERMOCOUPLE 40 


Ba 


-0.2480954 


THERMOCOUPLE 41 


Ba 


•0.3487988 


THERMOCOUPLE 42 


Ba 


•0.3038084 


THERMOCOUPLE 43 


Ba 


-0.3740253 


THERMOCOUPLE 44 


B» 


-0.94921B7 


THERMOCOUPLE 45 


Ba 


•1.0715818 


THERMOCOUPLE 44 


Ba 


-0.9038085 


THERMOCOUPLE 47 


Ba 


-0.8735839 


THERMOCOUPLE 48 


B* 


-0.B721191 


THERMOCOUPLE 49 


B- 


-0.8971479 


THERMOCOUPLE 50 


B- 


-0.4405957 


THERMOCOUPLE 51 


B« 


-0.8805175 


THERMOCOUPLE 52 


B- 


-0.40B442S 


THERMOCOUPLE 53 


B- 


2.2940912 


THERMOCOUPLE 54 


Ba 


-0.4140444 


THERMOCOUPLE 55 


Ba 


-0.4994428 


THERMOCOUPLE 54 


8» 


-0.244SB00 


THERMOCOUPLE 57 


B« 


-0.3954054 


THERMOCOUPLE 58 


Ba 


0.3308593 


THERMOCOUPLE 59 


Ba 


-0.4279294 


THERMOCOUPLE 60 


B- 


-0.0544453 


THERMOCOUPLE 41 


B- 


-0.2947773 


THERMOCOUPLE 42 


Ba 


-0.1212B90 


THERMOCOUPLE 63 


B* 


-0.7454101 


THERMOCOUPLE 64 


B- 


-0.4801757 


THERMOCOUPLE 65 


Ba 


-0.420B984 


THERMOCOUPLE 64 


Ba 


-0.1928710 


THERMOCOUPLE 67 


Ba 


-0.420B9B4 


THERMOCOUPLE 68 



Ma 


1.0028639 


Ba 


-0.5545898 


Ma 


1.0017774 


B» 


-0.4930175 


Ma 


1.0019627 


B- 


-0.4770995 


Ma 


1.0019207 


Ba 


-0.5316894 


Ma 


1.0026245 


Ba 


-0.5033203 


Ma 


1.0024214 


Ba 


-0.7976074 


Ma 


1.0017195 


Ba 


-0.4256835 


Ma 


1.0013046 


Ba 


-0.7514401 


Ma 


1.0019274 


B* 


-0.7330077 


Ma 


1.0026417 


Ba 


-0.9266112 


Ma 


0.9997830 


Ba 


-0.2324440 


Ma 


1.0019827 


B« 


-0.8003417 


Ma 


1.0017757 


Ba 


-0.4931152 


Ma 


1.0030603 


Ba 


-0.9011718 


Ma 


1.0023384 


Ba 


-0.7045917 


Ma 


1.0024405 


Ba 


-0.4196777 


Ma 


1.0010164 


Ba 


-0.3534133 


Ma 


1.0013094 


Ba 


-0.5520507 


Ma 


1.0013409 


Ba 


-0.4940937 


Ma 


1.0013094 


B- 


-0.5520507 


Ma 


1.0019588 


Ba 


-0.576B066 


Ma 


1.0011234 


Ba 


-0.S67B222 


Ma 


1.00171B5 


Ba 


-0.5253417 


Ma 


1 . 0232902 


Ba 


-0.B309081 


Ma 


1 . 00262B3 


Ba 


-0.7043456 


Ma 


1.0027475 


B- 


-0.6196288 


Ma 


1.0014108 


B- 


-0.4912597 


M» 


1.0025404 


B- 


-0.7730468 


Ma 


1.0019932 


B- 


-0.5812500 


Ma 


1.0020733 


Ba 


-0.471B749 


Ma 


1.0018148 


B» 


-0.59804 68 


Ma 


1.0018148 


B» 


-0 . 59B0468 


Ma 


1 . 001B339 


Ba 


-0.7208007 


Ma 


1.0009155 


Ba 


-0.5206543 



Two plots 
presented 
1 i n e a r i ty 



of T actual vs Tmeasures for 
in Figure A.l and Figure A. 2 



two thermocouples are 
to demonstrate their 
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APPENDIX B 



TABULATED DATA 

Heated Cylinder in the First Row 



H = 1.0 1.5m/s 3.0m/s 

Downstream Row No. Theta Theta 

1 1.0 1.0 

2 .468 .524 

3 .336 .380 

4 .251 .311 

5 .181 .211 

Heated cylinder in the first row. 

H = 4 . 6 1.5m/s 3.0m/s 

Downstream Row No. Theta Theta 

1 1.0 1.0 

2 .497 .484 

3 .365 .380 

4 .276 .277 

5 .197 .185 

Heated cylinder in the second row. 

H = 2 . 3 1.5m/s 3.0m/s 

Downstream Row No. Theta Theta 

1 1.0 1.0 

2 .501 .514 

3 .324 .354 

4 .250 .269 
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TABULATED DATA 



s,o(°C) 


T (° C ) 


Q ( w ) 


h ( w / m 2 ° c ) 


Re H 


N u 


U (m 


98.5 


20.3 


3.03 


56.61 


394.2 


14.09 


1.5 


98.9 


20.3 


3.05 


56.71 


394.2 


14.11 


1.5 


99.6 


20.5 


3.28 


61 .16 


393.9 


15.21 


i .5 


97.6 


20.0 


3.18 


60.34 


394.5 


15.02 


1.5 


96.3 


20.1 


3.74 


73.87 


394.4 


18.39 


1.5 


101.2 


20.8 


4.59 


87.27 


787.2 


21.69 


3.0 


89.4 


20.7 


4.17 


93.26 


787.4 


23.19 


3.0 


78.9 


20.5 


3.73 


98.52 


787.9 


24.51 


3.0 


96.8 


20.2 


4.74 


95.26 


787.6 


23.71 


3.0 


86.2 


20.3 


4.57 


107.69 


788.3 


26.80 


3.0 


111.0 


20.0 


4.12 


45.77 


907.4 


26.21 


1.5 


105.7 


20.1 


4.04 


• 47.93 


907.1 


27.44 


1.5 


102.8 


20.1 


4.13 


51.03 


907.1 


29.22 


1.5 


93.7 


20.1 


3.78 


52.61 


907.1 


30.12 


1.5 


99.3 


20.2 


5.50 


73.19 


1813.7 


41.90 


3.0 


110.1 


20.1 


5.99 


69.86 


1814.2 


39.99 


3.0 


118.2 


20.2 


6.81 


73.20 


1813.7 


41.90 


3.0 


71 .8 


20.1 


3.91 


80.03 


1814.2 


45.82 


3.0 


96.7 


20.0 


3.27 


42.74 


1814.7 


48.96 


1.5 


97.3 


20.2 


3.27 


42.49 


1813.7 


48.65 


1.5 


104.8 


20.3 


3.71 


44.20 


1813.2 


50.60 


1.5 


100.7 


20.6 


3.00 


36.9 


1811.6 


42.22 


1.5 


91 .5 


20.6 


2.92 


41.14 


1811.6 


47.06 


1.5 


112.3 


19 .5 


5.60 


62.83 


3634.7 


72.04 


3.0 


99.6 


19.8 


5.02 


65.70 


3631.6 


75.28 


3.0 


95.6 


19.9 


4.80 


66.23 


3630.6 


75.87 


3.0 


91 .2 


20.1 


3.92 


56.91 


3628.5 


65.17 


3.0 


85.8 


20.1 


3.92 


61.93 


3628.5 


70.91 


3.0 
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